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1. 


INTRODUCTION 


This  report  describes  the  work  performed  at  Science 
Applications,  Inc.,  to  improve  the  air/ground  interface  cor- 
rection factors  foimd  in  Version  4 of  the  ATR  (Air  Transport 
of  Radiation)  code(^) . Previous  versions  o f ATR(2,3)  had 
rather  simplistic  models  based  cn  the  work  of  French^^/  which 
does  not  involve  results  based  on  transport  calculations . 
Versions  3 and  of  ATR  contain  parameterized  results 

based  on  calculations  by  Pace,  et  al , in  which  only  a 
fission  and  a 14  MeV  neutron  source  were  used  and  thus  the 
parameterization  in  ATR  is  only  weakly  source  energy  depen- 
dent. Concurrently  with  another  effort^^) , two-dimensional 
adjoint  radiation  transport  calculations  were  performed  with 
the  use  of  the  DOT  computer  code.  These  calculations  formed 
the  data  base  for  the  TWEDEE  (Tactical  Nuclear  Weapon  Edit  of 
Dose  in  Endoatmospheric  - EartE  Environment)  coJe . ThTs  code 
was  used  to  generate  revised  ATR  correction  factors  for  the 
air/ground  interface  corresponding  to  those  monoenergetic 
groups  used  in  ATR  for  neutrons,  secondary  gamma  rays  and 
promps  gamma  rays  for  various  burst  heights  and  horizontal 
ranges.  It  should  be  noted  that  in  addition  to  improved  air/ 
ground  correction  factors  for  neutron  sources , this  effort 
represents  the  first  time  that  correction  factors  in  ATR  for 
gamma  ray  sources  (prompt  and  delayed)  are  based  on  transport 
calculations . 

The  air  over  ground  radiation  environments  computed  by 
Pace  and  the  TWEDEE  data  base  are  the  two  most  extensive  sets 
of  data  available  which  have  utilized  the  ENDF-IV  cross  sec- 
tion evaluations,  and  which  have  received  acceptance  and 
utilization  by  the  industry.  Besides  differences  in  source 
and  receiver  parameter  ranges , the  biggest  difference  in  the 
two  data  sets  is  the  type  of  ground  used  in  the  computations. 
The  Pace  calculations  utilized  what  is  referred  to  as  Nevada 
Test  Site  (NTS)  soil  and  the  TWEDEE  data  base  was  generated 
with  a dry  area  averaged  West  German  soil.  The  elemental 
constituents  of  these  two  soils  are  provided  in  Table  I.  The 
different  concentrations  of  hydrogen  in  the  two  soils  will 
have  the  largest  effect  on  the  radiation  transport.  In  a re- 
cent study  of  the  sensitivity  of  various  effects  to  the  radi- 
ation dose  from  nuclear  weapons  it  was  shown  that  a reason- 
able variation  in  water  content  could  result  in  differences 
in  the  dose  by  as  much  as  50%(^) . Thus  in  addition  to  the 
improved  capability  to  treat  source  spectrum  variation  and  to 
correct  for  the  air-ground  perturbation  for  photon  sources 
using  transport  data,  the  new  air  ground  correction  factors 
also  provide  another  soil  hydrogen  content. 
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Table  I.  Soil  Compositions  Used  by  Pace 
and  the  TWEDEE(^)  Data  Base 


Atomic  Concentration  (atoms/b-cm) 


Element 


(ATR  4) 


9.7656C-3)* 
3 ,4790(-2) 


4.8828C-3) 

1.1597(-2) 


TWEDEE  ^ ‘ ^ 

West  German  Soil 
(ATR  4.1) 

1, 

■10(-3) 

1, 

,02(-3) 

1, 

,90(-2) 

9, 

.9K-5) 

1, 

,61(-4) 

4. 

,36(-3) 

1 6, 

,83(-3) 

3. 

08(-4) 

9. 

,56(-4) 

8. 

,82(-6) 

1. 

■15(-6) 

1. 

■75(-4) 

Read  as  9 . 7656  x 10 


Table  II.  Source-Detector  Parameter  Ranges  for  the 
Air  Over  Ground  Data  Bases 


Pace 

Calculations 

T17EDEE  Data  Base 

Source  Energy 

2 

15  - 18 

Source  Height 

4 

14  (see  Table  IV) 

Detector  Response 

2 

Detector  Height 

4 selected 

1 

As  mentioned  previously,  the  two  sets  of  air  over 
ground  radiation  environment  data  bases  have  different  para- 
meter sets  for  the  source  and  target.  The  number  of  vari- 
ables for  the  source  energy  and  height,  and  target  response 
and  height  are  given  in  Table  II.  The  TWEDEE  data  emphasizes 
the  source  variation  but  has  only  one  target  height.  The 
Pace  data  has  more  information  on  the  target  response.  The 
work  reported  here,  concerns  the  efforts  to  combine  these  two 
different  and  somewhat  incomplete  data  bases  into  a form  suit- 
able to  provide  correction  factors  for  the  ground  perturba- 
tion for  all  source  spectra,  source  heights,  target  responses 
and  target  heights  encountered  in  the  analysis  of  weapon 
radiation  effects . 

The  correction  factors  were  parameterized  as  a func- 
tion of  horizontal  range  for  each  energy  group  and  each 
source  height,  and  the  resulting  coefficients  are  assembled 
into  Version  4 of  ATR  resulting  in  Version  4.1  of  ATR.  Ver- 
sion 4 remains  unaltered. 

For  the  delayed  photons  there  is  no  source  energy 
dependence  since  the  source  structure  is  inherently  build  in, 
thus  the  correction  factors  only  depend  on  the  geometry  con- 
figuration . 

This  report  describes  the  improved  air/ground  data 
base,  its  parameterization  and  subsequent  results  of  ATR 
Vers  ion  4.1. 


2. 


BACKGROUTJD 


The  basic  quantity  constituting  the  data  base  for  ATR 
is  the  homogeneous  infinite  air  fluence  $g (Eg , Eq , o , u ) which 
is  a function  of  source  energy  Eg,  detector  energy  Eg,  range 
c in  units  of  gm/cm"-  and  angle  u at  the  detector.  3y  proper 
scaling  techniques  the  air  density  variations  can  be  accounted 
for  at  arbitrary  geometry  configurations,  but  other  models  are 
necessary  to  correct  for  the  air-ground  interface  when  the 
source  and/or  the  target  are  within  about  one  kilometer  of  the 
ground. 

(0  3-) 

The  initial  versions  of  ATR  ’ used  variations  of  the 
first-last  collision  model  of  French^S).  From  this  model 
factors  F(H^)  and  F(K’g)  are  derived  for  source  height  He  and 
target  heignt  Hy,  respectively,  which  are  treated  as  indepen- 
dent perturbation  factors  to  correct  the  infinite  air  results 
according  to  the  following  equation: 

t(E2,Eg,p,H3,H.j.)  = OQ(E2,Eg,p)F(H2)F(S.j)  (1) 

where  the  correction  factors  introduced  two  independent  vari- 
ables , Hq-  and  the  perturbation  in  the  angular  dependence 
was  not  considered.  The  first-last  collision  model  estimates 
F(Hg)  by  calculating  the  effective  fraction  of  first  colli- 
sions about  the  source  as  a function  of  source  altitude,  Hg . 
The  function  F(Ht)  is  determined  by  calculating  the  number  of 
last  collisions  in  the  vicinity  of  the  target. 

In  addition  to  this  model,  for  neutrons  and  secondary 
gamma  rays,  the  correction  factors  were  modified  by  Straker's 
two-dimensional  transport  calculations ( 8)  of  dose  versus 
distance.  For  prompt  gamma  rays  only  a slightly  modified 
first-last  collision  model  exists  for  all  ATR  versions  . 


Pace  at  Oak  Ridge  National  Laboratory  performed  two- 
dimensional  radiation  transport  calculations  in  air-over- 
ground geometries (6)  with  the  DOT- III  code  for  14  MeV  neutron 
source  and  a weaponized  fission  source.  The  data  for  neutrons 
and  secondary  gamma  rays  was  reported  for  four  different 
source  heights:  1,  50,  100,  300  m,  a set  of  target  heights  to 
approximately  1300  m and  a set  of  horizontal  ranges  to  about 
1500  m.  This  data  was  obtained  and  the  following  correction 
ratio  was  parameterized  within  ATR  Version  4 : . ) 
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(2) 


i 


i 


where  R is  the  horizontal  range,  Hg  is  the  source  height,  H-p 
is  the  target  height,  Dq(R,H5  ,H'p)  represents  the  Pace  Hender- 
son tissue  response  and  D^CR.Hg.H-p)  represents  the  corre- 
sponding ATR  infinite  air  data.  The  data  is  source  energy 
dependent  only  to  the  extent  that  two  sources , 14  MeV  and 
fission,  were  used  when  the  data  base  was  generated.  The  cor- 
rection factors  are  used  to  multiply  the  ATR  infinite  air 
results  to  obtain  the  air-over-ground  effects.  For  source 
energies  below  5 MeV  the  fission  source  correction  factors 
are  used,  and  for  source  energies  above  5 MeV  the  14  MeV 
source  correction  factors  are  used.Cl-'^) 


Thus,  ATR  Version  3 and  Version  4 generate  fluence  and 
dose  quantities  described  by  the  following: 


KE2,E^,R,Hg,H.^)  = C(R,H5,H.^)  (3) 

Of  course,  the  horizontal  range  R and  slant  range  t are  re- 
lated because  of  the  geometry  configuration.  Also,  inherent 
in  the  above  equation  is  that  the  source  energy  dependence  is 
limited  to  the  5 MeV  cross  over  selection  constraint  previ- 
ously described. 
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3.  SOURCE  ENERGY  DEPENDENCE  OF  AIR-OVER-GROUND  INTERFACE 


A further  improvement  of  ATR's  air-over-ground  inter- 
face correction  factors  is  obtained  by  utilizing,  a computer 
code  called  Tl'J'EDEE  (Tactical  Nuclear  Edit  of  Dose  in  Endo- 
atmospheric-Earth  Environment)  which  was  deveToped  at  Science 
Applications,  Inc..  (7)  This  code  folds  source  energy  im,por- 
tances  and  sources  for  prompt  and  delayed  dose  for  detectors 
located  at  1.75  m from  the  ground  using  a data  base  generated 
with  the  DOT  code  in  ajoint  mode.  The  soil  composition  used 
in  the  Face^^)  forward-mode  DOT  calculations  and  the  soil 
composition  used  for  the  D0T/TUEDEE(^)  adjoint  transport/data 
base  generation  are  listed  in  Table  I. 


Using  the  Tl'JEDEE  code,  it  is  possible  to  obtain  the 
following  dose  values:  D'pCEg , R,Hg  ,Ht  = 1.75  m)  for  neutrons, 
secondary  gamma  rays  and  prompt  gamma  rays  for  each  mono- 
energetic  source  band  corresponding  to  the  ATR  structure  for 
a set  of  horizontal  ranges  to  1300  m and  source  heights  to 
1100  m.  Corresponding  infinite  air  dose  values  are  calculated 
by  ATR  as  well,  and  the  following  correction  ratio  is  formed: 


C,p(Eg  , R,Hg  ,H,p  = 1.75  m) 


D.j,(Eg  .R.Hg  ,H~  = 1.7  5 m) 
D,  (Ec  ,K..„  = 1.75  m) 

Ad  o i 


(4) 


In  order  to  account  for  the  target  height  variation,  the  cor- 
rection factor  is  multiplied  by  the  following  ratio; 


P(R,H2,H.^) 


C ( R,  , H.J,) 

C(R,Hg,H,p  = 1.7  5 ml 


(5) 


For  neutrons  and  secondary  gamma  rays  this  perturbation  factor 
is  available  from  the  Pace  data.(o)  Note  that  the  correction 
for  target  height  is  based  on  data  from  calculations  performed 
with  a different  hydrogen  concentration  in  the  soil  than  the 
data  base  used  for  Eq . (4).  Although  this  procedure  is  less 
than  ideal,  it  is  felt  to  be  more  reliable  than  other  schemes 
such  as  a last  collision  model.  Again,  the  attempt  is  to  use 
the  information  available  to  produce  the  most  accurate  model. 
Figures  1 through  d represent  exam;ples  of  this  perturbation 
factor  for  neutrons  and  secondary  gamma  r.ays  , one  figure  for 
each  source.  There  perturbation  factors  are  used  to  correct 
for  the  target  height  variation  for  neutrons  and  secondary 
gamma  rays.  For  the  prompt  and  the  delayed  gamrr.a  rays  the 
first-last  collision  mode'l*'^^  is  used  to  correct  for  target 
height  variations. 


I 

I 


I 

I 
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Pe  rtiurbat  ion  raiio  factors  for 
a fission  neutron  source  at  100 
target  heights . 


4. 


GENERATION  OF  THE  NEW  ATR  AIR/ GROUND  DATA  BASE 


Dose  responses  for  neutrons,  secondary  gamma  rays  and 
prompt  gamma  rays  are  generated  from  the  TWEDEE  code,  the 
scattered  components  of  each  are  selected  and  the  correction 
ratio,  Cij,  (Equation  4),  is  generated. 

Table  III  contains  the  selected  horizontal  ranges  and 
fable  IV  contains  the  selected  source  heights  at  which  the 
data  is  generated.  The  choice  of  these  parameters  is  dic- 
tated by  the  possible  range  of  the  TWEDEE  data  as  well  as  the 
areas  of  anticipated  structure  of  the  correction  factors. 

Tlie  energy  group  structure  in  the  TWEDEE  code  is  dif- 
ferent from  that  of  the  ATR  code,  thus  some  regrouping  is 
necessary  and  some  of  the  ATR  energy  groups  are  split  into 
several  of  the  TWEDEE  energy  groups . The  total  span  of  the 
energy  structure  of  TWEDEE  includes  the  total  span  of  ATR. 

The  correction  factor  data  base  was  generated  by  form- 
ing the  ratio  of  the  TI'JEDEE  dose  to  the  ATR  4 infinite  air, 
uncorrected  dose  for  its  given  source- target  configuration. 

The  neutron  dose  response  used  for  the  ratio  was  the  tissue 
kem*a  and  photon  response  was  the  Henderson  tissue  dose. 

Figure  5 shows  a representative  sample  of  the  new  neu- 
tron correction  factors  for  different  source  energies.  Due 
to  certain  convergence  limitations  of  the  DOT  code  some  of 
the  neutron  data  corresponding  to  low  energies  is  not  avail- 
able from  T^NEDEE.  This  limitation  is  not  very  serious  for 
most  applications  since  it  occurs  for  source  heights  of  500 
meters  and  above,  for  source  energies  of  3.35  x 10~3  MeV  and 
below.  A sensitivity  analysis  shows  that  for  a typical  thermo- 
nuclear source  those  low  energy  sources  contribute  less  than 
1.5%  to  the  total.  For  a typical  fission  source  those  low 
energy  sources  do  not  contribute  at  all. 

Figure  6 shows  a representative  sample  of  the  new  cor- 
rection factors  for  neutron  generated  secondary  gamma  rays . 
Experience  has  shown  that  secondary  gamma  ray  dose  variation 
is  not  as  smooth  as  the  neutrons  or  prompt  gamma  rays  varia- 
tion, thus  it  takes  more  coefficients  to  represent  each  curve 
as  it  will  be  detailed  in  the  next  section.  The  data  for 
secondary  gamma  rays  is  everywhere  dense  in  that  TWEDEE  data 
exists  for  all  source  heights  and  source  energies. 

Figure  7 shows  examples  of  the  new  correction  factors 
for  prompt  gamma  rays.  As  expected,  the  curves  exhibit  a 
generally  friendly  behavior  and  each  may  be  parameterized  by 
fewer  coefficients  than  neutrons  or  secondary  gamma  rays. 
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Figure  6.  Secondary  gamma  ray  correction  factors  as  a 
function  of  horizontal  range  for  various 
neutron  monoenergetic  sources  and  a source 
height  of  1 meter  and  a target  height  of 
1 . 75  meters . 
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The  TWEDEE  group  structure  below  0.2  MeV  was  found  to 
be  too  broad  to  give  accurate  results  for  source  photons  be- 
low this  energy.  As  detailed  in  Ref.  2,  it  is  necessary  to 
partition  the  energy  range  below  0.2  MeV  into  many  groups 
(on  the  order  of  30)  to  get  reliable  data.  Comparisons  be- 
tween fine  group  and  broad  group  results  showed  discrepancies 
as  large  as  a factor  of  two  in  this  regime. (2)  However,  a 
sensitivity  analysis  was  performed  which  showed  that  the  con- 
tribution to  the  total  dose  from  source  photons  of  energy 
less  than  0.2  MeV  is  on  the  order  of  1 . 57o  for  reasonable 
source  spectra  at  ranges  of  interest.  It  was  decided  to  use 
the  existing  scheme  of  the  first-last  collision  model(l)  to 
correct  the  ATR  results  for  source  energies  below  0.2  MeV. 

The  fission  product  gamma  ray  correction  factors  were 
generated  by  using  the  source  values and  integrating  over 
the  source  spectrum,  thus  the  correction  factors  are  only  geo- 
metry dependent,  not  energy  dependent.  In  order  to  alleviate 
the  disparity  of  the  lower  energy  groups  below  0.2  MeV,  the 
source  values  corresponding  to  those  energies  are  not  used  in 
the  computation  of  the  correction  factors . 
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5.  PARAMETRIZATION  OF  AIR/GROUND  CORRECTION  FACTORS 


Three  subroutines  (AGNEUT,  AGGAMM,  AGSGAii)  have  been 
written  and  tested  that  provide  the  parametrized  air/ground 
correction  factors  for  neutrons , prompt  gammas  and  secondary 
gammas.  These  routines  use  coefficients  generated  by  a piece- 
wise  least  squares  fitting  of  the  appropriate  data  to  alge- 
braic functions  of  various  modified  polynominal  forms.  The 
parameters  of  the  fits  are  heuristically  selected  so  as  to 
obtain  the  minimum  error  and  least  number  of  coefficients. 

The  coefficients  so  generated  are  arranged  by  source  energy 
groups  for  ease  of  use  and  incorporation  into  the  ATR  code. 

For  the  fission  product  gamma  rays  a set  of  coefficients  are 
generated  for  each  source  height  and  fit  as  a function  of 
horizontal  range. 


The  basic  operation  of  the  functions  is  as  follows. 
First,  the  source  height  desired  is  delimited  between  the 
points  of  the  data  base  after  which  the  desired  horizontal 
range  is  isolated.  These  values  are  then  used  to  locate  a 
pointer  word  which  contains  information  about  the  location 
and  number  of  coefficients  and  type  of  function  to  be  used  in 
computing  the  air/ground  ratio.  The  ratio  is  then  computed 
for  adjacent  source  heights  and  a linear  interpolation  between 
source  height  values  is  used  to  return  the  final  value  of  the 
desired  correction  factor.  When  the  horizontal  range  is  less 
than  three  meters,  an  extrapolation  is  used  to  assign  an 
appropriate  air/ground  ratio  by  using  the  results  at  three 
meters  and  at  five  meters . 


The  modified  algebraic  functions  used  to  parametrize 
the  original  data  base  assume  the  following  forms,  where 
is  the  ratio  for  neutrons  and  prompt  gamma  rays  is  the  log 
of  the  ratio  for  secondary  gamma  rays . 


Function  Tvpe 

_0  ; 

^ ^ 3 

a^  X + a^  X 

2 ^ 1 ^ 

+ a2  X + a^  X + ag  = 

Ct 

(6) 

Function  Type 

J.  ; 

3 2 

a/  X + a,  X 
u p 

+ a^  X + a^/x  + ag  = 

Ct 

(7) 

Function  Type 

2 ; 

a^  /x'^  + a^/ x^ 

•4  J 

9 

+ a.,/x^  + a^/x  + ag  = 

Ct 

(8) 
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Function  Type  3 ; 


3 2 - 

X + X + a^  X + a^^vx  + a^  = 


(9) 


In  addition,  an  extrapolation  function  is  used  when  the  hori- 
zontal range  is  less  than  three  meters . If  the  horizontal 
range  is  larger  than  the  span  of  the  data  base,  then  the  value 
corresponding  to  the  last  horizontal  range  in  the  data  base  is 
used . 


The  pointer  word  which  contains  information  regarding 
the  location  and  number  of  coefficients  and  the  type  of  func- 
tion is  composed  as  follows , numbering  from  the  right  with  the 
rightmost  binary;  bit  numbered  one. 


Source  Height  Index 
Horizontal  Range  Index  (start) 
Horizontal  Range  Index  (stop) 
Number  of  Coefficients 
Type  of  Function 

Starting  Location  of  Coefficients 
Within  Coefficient  Array 


■ woras  are  grouped  according  to  source 
energy  group.  Within  each  group  there  is  a pointer  word  for 
each  piecewise  fit.  For  example,  assume  that  the  third  through 
eighth  points  of  source  energy  group  six,  source  height  ten 
were  fitted  to  the  function  of  the  form  a^^x  + an  = C ^ . The 
pointer  word  containing  this  information  would  be  located  in 
index  group  six  (for  source  energy  group  six)  and  appear  in 
binary  as  follows; 


Bits 

1- 

/ 

^4 

Bits 

5- 

8 

Bits 

9- 

12 

Bits 

13- 

15 

Bits 

16- 

17 

Bits 

18- 

28 

Thest 

; DO 

in 

Location  of  Coef s . 

00  1 

010 

1000 

0011 

1010 

28  18  16  13  9 5 1 


For  machine  to  machine  compatibility  the  pointer  words  have 
been  stored  as  integers  using  "DATA"  statements,  thereby 
obscuring  the  structure.  However,  a simple  integer  to  binary 
conversion  will  reveal  the  format. 

The  coefficient  array  consists  of  all  coefficients 
generated  by  the  fitting  process  arranged  in  the  order  of 
source  energy  group,  source  height,  horizontal  range  and  is 
ever}A>7here  dense.  The  starting  location  of  the  coefficients 
within  the  array  for  a particular  piecewise  fit  is  stored  in 
the  corresponding  pointer  word  as  previously  described. 
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This  process  resulted  in  a total  of  3365  coefficients , 
1073  to  generate  the  neutron  data  base,  1249  to  generate  the 
prompt  gamma  data  base  and  1043  to  generate  the  secondary/ ^ 
gamma  data  base.  The  original  data  base  contained  approxi- 
mately 9800  points.  Therefore,  the  parametrization  resulted 
in  a reduction  of  approximately  6435  data  points  or  65  per 
cent. 
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6.  ATR  VERSION  4.1 


j 


The  routines  described  in  the  previous  section  were  ; 

incorporated  into  ATR  Version  4(^)  and  the  new  version  will 
be  called  4.1.  Here  the  results  are  presented  for  each 
"particle”  type  at  several  geometry  configurations . There 

are  four  figures  for  each  radiation  type  for  a 1 KT  burst,  i 

each  figure  containing  three  curves : the  first  curve  is  with-  - 

out  air/ground  correction  factors  (scaled  infinite  air  re-  j 

suits),  the  second  curve  shows  the  results  using  the  cor- 
rection factors  present  in  ATR  Version  4(4)  and  the  third 
curve  shows  the  results  due  to  the  use  of  the  new  energy 
dependent  correction  factors  (.ATR  Version  4.1). 

The  four  geometry  configurations  are  selected  in  a 
manner  that  shows  the  effects  of  both  the  source  and  the  tar- 
get close  to  and  away  from  the  ground.  The  first  figure  in 
each  sequence  has  both  the  source  and  the  target  at  one  meter 
above  the  ground,  the  second  one  has  the  source  at  one  meter 
and  the  target  at  800  meters,  the  third  one  has  the  source  at 
800  meters  and  the  target  at  one  meter,  and  finally,  the 
fourth  one  has  both  the  source  and  the  target  at  800  meters. 

Figures  3 through  11  show  the  results  for  neutrons; 
the  source  used  is  a typical  thermonuclear  spectr'cm  with  a 
normalization  of  1.9  x 1023  neutrons/KT.  Figures  12  through 
15  show  the  results  for  prompt  gamma  rays,  the  source  used  is 
a typical  fission  spectrum  with  normalization  of  10^2  gammas/ 

KT . Figures  16  through  19  show  the  results  for  secondary 
gamma  rays  and  Figures  20  through  23  show  the  results  for 
fission  product  gamma  rays.  Note  that  the  largest  differences 
between  ATR  4 and  .ATR  4.1  results  occur  for  the  gamma  ray 
sources  (prom.pt  and  delayed)  . The  correction  factors  for 
gamma  ray  sources  for  ATR  4.1  are  based  on  transport  data 
while  the  ATR  4 corrections  were  based  on  an  unverified  model. 

There  is  no  visible  impact  to  the  user  as  far  as  the 
ATR  Version  4.1  command  structure  is  concerned,  only  some  of 
the  internal  mechanism  is  different.  The  Appendix  contains  a 
composite  sample  nroblem  that  illustrates  the  use  of  ATR 
Version  4.1. 

The  ATR  4.1  version  has  the  follow  features  which  are 
not  present  in  ATR  4: 

• .Air  ground  correction  factors  which  explicitly 
treat  variations  in  the  neutron  and  gamma  ray 
source  energy  spectra 
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Figure  12.  Prompt  gamma  ray  Henderson  tissue  dose  for  a 1 KT 
yield,  source  height  of  1 meter,  target  height  of 
1 meter  and  various  air/ground  correction  factors 
for  a typical  fission  source 
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gure  13.  Prompt  gamma  ray  Henderson  tissue 
yield,  source  height  of  1 meter, 
800  meters  and  various  air/ground 
tors  for  a tvnical  fission  source 
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Figure  15.  Prompt  gamma  ray  Hender.son  tissue  dose  for  a 1 KT 
yield,  source  height  of  800  meters,  target  height 
of  800  meters  and  various  air/ground  correction 
factors  for  a typical  fission  source. 


SKCONDARY  GAMMA  RAY  HRNDRRSON  TISSUE  DOSE  (RADS) 


Figure  18.  Secondary  gamma  ray  Henderson  tissue  dose  for 
1 KT  yield,  source  height  of  800  meters,  targe 
height  of  1 meter  and  various  air/ground  correc- 
tion factors  for  a typical  thermonuclear  neutron 
source . 
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Figure  21.  Fission  product  gamma  ray  Henderson  tissue  dose 
for  a 1 KT  yield,  source  height  of  1 meter,  tar 
height  of  800  meters  and  various  air/ground  cor 
rection  factors  for  a delayed  gamma  rav  source. 
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Fission  product  gamma  ray  Henderson  tissue  dose 
for  a 1 KT  yield,  source  height  of  800  meters, 
target  height  of  300  meters  and  various  air/ground 
correction  factors  for  a dealved  gamjna  rav  source . 
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• Air  ground  correction  factors  generated  from 
transport  calculations  in  which  the  ground 
composition  is  more  realistic  for  tactical 
scenarios  than  the  NTS  composition 

• Air  ground  correction  factors  for  gamma  ray 
sources  which  are  based  on  detailed  transport 
calculations  rather  than  a mathematical  model. 

These  additional  features  will  permit  a more  accurate  analysis 
of  the  radiation  environments  in  tactical  nuclear  weapon 
scenarios . 


It  should  be  noted  that  there  is  no  set  of  correction 
factors  which  could  be  incorporated  into  ATR  which  would  oe 
appropriate  for  all  possible  analysis  of  weapon  effects  in 
air  over  ground  configurations.  The  number  of  combinations 
of  source  height,  source  energy,  target  height,  target  energy, 
target  angular  distribution,  and  ground  composition  is  too 
large  to  even  hope  to  generate  a data  base  or  to  incorporate 
such  a data  base  into  ATR.  There  are,  however,  three  remaining 
portions  of  the  air  over  ground  radiation  environments  which 
could  stand  reasonable  improvement  with  models  based  on  calcu- 
lations. These  are: 

• Variation  of  the  dose  as  a function  of  hydrogen 
content  in  the  soil 


• Variation  of  the  dose  as  a function  of  target 
height 

• Estimation  of  the  "upward”  and  "downward"  component 
with  respect  to  the  ground  of  the  radiation 
environment . 
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APPENDIX  A 


SAMPLE  PROBLEM 


This  sample  problem  is  a composite  of  four  different 
radiation  types:  neutrons,  prompt,  secondary  and  fission  pro- 
duct gamma  rays.  The  source  selection  for  neutrons  and 
secondary  gamma  rays  is  a typical  thermonuclear  spectrum  with 
a normalization  of  1.9  x 10^3  neutrons,  for  prompt  gamma  rays 
it  is  a typical  fission  spectrum  with  a normalization  of  10^3 
gamma  rays.  The  fission  product  yield  is  specified  as  1 KT . 

The  geometry  configuration  is  given  as  the  source  and 
target  one  meter  off  the  ground  and  the  horizontal  range 
varies  from  10  meters  to  1400  meters. 

The  output  request  is  simply  all  of  the  dose  responses 
for  all  of  th®  radiation  components. 

The  output  contains  the  details  of  the  source  con- 
figurations followed  by  the  output  of  the  dose  response  for 
all  of  the  requested  components. 
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A.l  SAMPLE  PROBLEM  INPUT  COMMANDS 


*TiTLE  SAmPi  F ppopih'vi  F.-)R  ATW  VFPSt''''''  . 1 

*N-SOU»CF (?) 

*N-(^r)NM  1 ,q*.?3 

♦G-SOUPCF  f 1 ) 

*G-Nf)PM  i+?3 
*FP-Y  t 
1 

*MT  1 

*RH  10  50  too  200(^00)1000 

*r)fJSE/'^/ 

*onsE/G/ 

*i;nsF/MG/ 

*DnsE/Fp/ 

♦EXC 

*FIM 
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SAMPLE  PROBLEM  OUTPUT 


A. 2 


VE9SI0M  0,1  np  ATPMxr;  - phl  «iR/r.wnijNn 

4TR  PwnnLEM  MIIHRFR  1 SARPLF  prohle**  EHb  jtR  vpRStnR  a.l 


reijtrhn  sniiPCF  r*jrFRN*L  thfrhoukClfar 
normal  1 inn=l  .TOOt  NfilTRON  /KT,  Y I p LOs  1 . n n OF  ♦ rt  0 KT 
TOTAL  OMTPIITil  NfUTRON 


pNpRljYtMpv  1 

1 . nTE-nS-^.RnF-OR 

2, Rnp-n.;-,^n,P-n« 
1 .0  1 F -OU'5  .RTE-O'I 
5.fl3E-oa-3, T5E-0^ 
3. 3SE-0  5"J • 1 1 ^"0 1 
1 .t |E-ni-S.S0E-Pl 
'5,'5nE-0|*i,]|E*0ft 
1 .1 lF*no-|  .«3E*on 
1 .S3E*on-2.35E+on 


SOURCE  SPFCTRUM 

N N/MpV  pNpRr.Y(Mpv  ) 

o.finE-ni  p.noE-01  ?.5SF*no-?.'J'>F*on 


3 . roe 

« .ShF » ->| 

1 .R«E»?? 
! 

t . I SEfja 
S.3?F*?l 


S.?«E  + 2/i 
P . 06E  + ?u 
B.3''F*2a 
M.auF  +23 
«.aiE  + ;>? 
2.««E*?? 
1 .ftaF+?? 
1 .P?F ♦?? 


j.ahF«-nn-5,P|P+0P 
3.PiE*PO«a.n7FYpp 
0.07E»00-0,R7F+PP 
<i.P7F*on-h,3''F  + pp 
6.36E+P0-R, lOF*Pp 
R . 1 RP  + PO-  1 , OPF  + IT  I 
1 ,0pF*PI-1 .??F+01 
1 .2?F  AO  1-1  .SPFaO I 


N 

R.SpFaJp 

T.BiEa?! 

n.P'lFf^l 

3.?3F+?I 

3.R2Fa?i 

?.70Fa?i 

J.^bEa?, 

0 , RrE ♦? 1 
1 


fJ/Uf  V 
R.R'lF  »?1 
R.SRF  Aj>l 
P.RRF  »J1 
3.SOF  A?1 
?.'I0F  *?1 
1 .53E  a2| 
1 .RBE*?! 
?.?lE  »?1 

'I  . 7RF  A^l 


ATR  PRORlEh  nUmRFR 


sample  PRt)R|  Em  F|\r  atR  VERSION  a.l 


gamma  Source  iNlERfAL  FISSTOn 

NOP"AL  I 7A  T TONS  1 . nopp  3 GAHM4  /K  T , Y I F I,  0 r I , P P 0 E A a P KT 
TOTAL  OUTPIII-l  .opof  a?3  gamma 

SOURCE  spectrum 


FNfRUT(mfV  1 

V 

pNFRGY(Mpv  ) 

G 

G/MfV 

? . PpE -07-5 . OOE -p ? 1 

l.oPE,.?! 

1 .03Fa23 

1 . 33E  aOO- 1 .hPF  aPO 

P. 1 RE  A?i 

1 ,R7Ea?7 

Opf -07- 1 , floE -p 1 1 

. 3SE  A71 

7. 7|Fa^P 

l.opFAnp-7.npFA/7n 

3.R3EA71 

1 . 1 RE  *77 

1 . OO'-P 1 -2 . OPF-P 1 F 

MPF+pi 

R,  1 PFa?7 

?#ooF^oO*’?.*'OF^nft 

3.  7a.E  A?1 

7.5|F  A?1 

2, OPE «P 1 -3, PPE-P  1 * 

..PTEaji 

A..R7E 

?.SPFapp-5,ppe*pp 

7.73'^»71 

p.a7E  Ap, 

3,  ooE-r,  1 -a , ooE-p  1 " 

'.R«E  a7I 

R.faPF  A7? 

3.OPFAPp>a,0OE+PP 

7. !7Fa?i 

7. 1 7f  *71 

«.OpF-oi-6.floE-01  1 

. 77E  a?7 

R.RaEA7? 

U.OftF  4(>'> 

7 , 'ARE  a7p 

7.0RFA7P 

<>.Opf"0I-'’,opF-p,  1 

.apEA?7 

7.piF*77 

S.PpE ♦PO-0,SpFAnp 

3.73EA7n 

7.  1 “iF  *7p 

R.POF-PI -1 .OpFaOP  I 

.PoE  »?? 

S.o?Fa7? 

O.SoEaPP-P.OPEaPP 

R.7oFAiP 

a.s3E air 

1 .00^*00*1  •33E  + PP  1 

.0  I’"*?? 

3.7'5F*7? 

R.ppF  Ypo'!  .OpEY-fli 

I ,SRFa |Q 

7.RpE*|R 

ATR  proMlFm  mUmmfr 

1 

SAmP|  F PRIlRi  F M F|1R 

ATR 

Tfuj  A , I 

FlSSln*)  PRonuCT  r r F|  Os  1 . OOF  apPr  T . 1 

fission  FUAriTtlNP  1 

1 .Anr^Ao 

ATR  PROH)  Fm  nUmRFR 


SAmPi  F PR<in|F«  Fur  atp  vfrrT|.n  i , I 


GRoiiNn  level 

mopIZ.  R'aiGF 
• SI  ANT  RANGF 


O.PPPK“, 
□ H-  p , p 1 Pa  M , 
«S-  P^PIPKM, 


p . pppr.M  /r  “• 
t ,??sr,u/rMA 

1 ,P?e;(.M/ru, 


•7#  P.PPPKET 
• 7 r P . p 3 3«  F T 
*7,  P.P33KFT 


, p.pppMiirs 

, n.pPNMliF*; 

, p.pp»,“n  ‘ F 
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A. 2 (continued) 


T 4RGF  T 

‘LT. 

HTs 

0 . 00  I K M r 

n. 

1??gm/Cm*»?,  n 

.oo'Skft  , 

O.OOlMltl-t^ 

S'lURCE 

4L7. 

HS  = 

0,001km, 

n. 

12?r,M/cH««?,  n 

^ nn  F T , 

o.ooiMjt.^^s 

• SL4»<T 

4NGLF 

AKJ. 

o.onnnFGBEES 

(cns=  i.noonni 

*r*Lr"L»TFn  f”""  niHFi' 


sFiiTPOM 

OOFF 

vs.  nn»rz.  DifiCF  rp*Ps  i 

nnSF 

t = 

'•lE'inFuSHN  TISSMF 

noFF 

? = 

SNVnFq.MEliFFI  D 

f'osF 

Tissue 

nnsE 

/I  5 

MI  r)-RM«NTriM 

^OSE 

s = 

CONF  BE  TE 

nn3F 

hs 

M» 

nnsE 

7 = 

nnFi-I0»JlZlNG  Si. 

DMSE 

«r 

rnNtzr‘'G  sn_TC''N 

HnPtZ,  PANHF 

nnSF  1 

OOSF  ? rniSF  S riOSF  a r>nsF  S OnSF  6 

nnsE  7 

nrisr 

1 . nnoF  »n 1 

M 

3."?Etn7 

b.3SF*07 

d.ZoF+07 

I .BEEfO? 

b.?BF  tOb 

B . 3bF  t Ob 

S , OpP  fOS 

Z.OSF  tOb 

5 . 0 n n ^ F fi  1 

M 

1 . SSE  Fob 

Z.U?E  + oT' 

1 .aBEFob 

b.hoE  F()G 

?.  I'iEfos 

Z . BbE  t riS 

1 ,73Ftoa 

b . b 0 F * 0 0 

1 . nfinE  tn? 

M 

3.Z5E..05 

p. 1 SF  *03 

S.hSFfoS 

1 ,S«FFflS 

S , 0 UE  t OB 

b.BSEtno 

n.p^F  *r^\ 

I , apf  1 0 u 

?.0OOE*f>c 

M 

h.07E*n<) 

1 .I?E+0S 

b . bBF  ♦ (Ml 

Z.SoF  *00 

1 oE  *03 

1 .zaEty/| 

s.  1 BEtoz 

Z . boF  to  3 

« . noOF+nZ 

M 

t."SF,n3 

1 .Z'’F  + 0(1 

7 .SbF  *0  3 

3.01Ff03 

B ,BsF+0? 

1 .3af  f03 

'7.3lF  tOl 

Z.30F  top 

b.OOOFtBz 

1 .?0EhT5 

Z.?bF»(l^ 

1 . 3F'’F03 

S.nSFFB;, 

I ,Tnr*np 

p.pp^ 

I . b TF  t 0 ( 

3.b  7F  to  , 

« .finoF  toz 

M 

Z.73F  *np 

S.Z3F  FO? 

3 . 03E  + 0? 

1 . 1 ?F  FO? 

3.«|Ftni 

il  . B 1 F to  1 

3.73Fton 

7 , ■'!  0 F ton 

1 .onOF  +n3 

M 

7.a5E*oi 

1 .aaF*o? 

«.?sF*m 

P , B«F  to  1 

1 .nZFtoi 

1 . 3 0 E 1 0 I 

1 .npE  ton 

1 . « 0 E t n n 

1 •Zor'E  *ni 

M 

?.3tE*ol 

U.SlF  *(11 

Z . SbF  F a 1 

B.  1 |F  trio 

3.1  |F  ton 

3.BBE  too 
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